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2 ' Abstract 

o 

We investigate the effect of a microscopic three-body force on the proton and 
neutron superfluidity in the 1 Sq channel in /3-stable neutron star matter. It is 
^ | found that the three-body force has only a small effect on the neutron Sq pairing 



gap, but it suppresses strongly the proton 1 So superfluidity in /J-stable neutron 
star matter. 
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Superfluidity plays an important role in understanding a number of astrophysical 
phenomena in neutron stars [1-10]. It is generally expected that the cooling processes 
via neutrino emission [5, 6, 7], the properties of rotating dynamics, the post-glitch 
timing observations [8, 9], the possible vertex pinning [10] of a neutron star are rather 
sensitive to the presence of neutron and proton superfluid phases as well as to their 
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pairing strength. For instance, since the paired nucleons do not contribute to thermal 
excitations, the proton superfluidity could suppress considerably the neutrino emis- 
sion processes, and consequently affects the neutrino cooling rate of a neutron star 
remarkably. 

Since the neutron and proton superfluidity properties in neutron stars are related 
only indirectly to the observations, reliable and precise theoretical predictions based on 
microscopic many-body approaches are highly desirable. In Refs. [11, 12] it is reported 
by using the Brueckner-Hartree-Fock ( BHF ) calculations based on purely two-body 
nucleon-nucleon ( NN ) interactions that the 1 So neutron superfluid can be formed 
only in the low-density region ( p B < 0.1 /m~ 3 ) with a maximal gap value of about 
2.8MeV at a Fermi momentum kp — 0.8fm _1 , while due to the small proton fraction 
in /^-stable matter the proton superfluid phase in the 1 So channel may extend to much 
higher baryon densities up to ~ 0.4fm -3 with a maximal pairing gap of about 
0.9MeV. In Ref. [13], the relativistic effects on the superfluidity in /^-stable neutron 
star matter have been investigated by using the Dirac-BHF ( DBHF ) approach and 
the one-boson-exchange iViV interaction ( BONN potential ) [14]. It was found that 
the relativistic effect on the proton superfluidity in the 1 S'o channel is very weak, while 
the magnitude of the neutron pairing gap in the 3 P2 channel is reduced strongly due 
to relativistic corrections to the single-particle ( s.p. ) energies. 

The pairing correlations in nuclear medium are essentially related to the underly- 
ing NN interaction and their magnitude is determined by the competition between the 
repulsive short-range and attractive long-range parts of the interaction. The solution 
of the gap equation is extremely sensitive to the medium modifications of both the 
bare iViV interaction and the s.p. energy. These medium effects on nuclear pairing 
are extremely difficult to study on a microscopic level and have motivated an exten- 
sive investigation by many authors [15, 16]. Three-body forces, which turn out to 
be crucial for reproducing the empirical saturation properties of nuclear matter in a 
non-relativistic microscopic approach [17, 18, 19], are expected to modify strongly the 
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in-vacuum AW interaction, especially the short-range part [17, 19, 20]. However, their 
effects on the superfluidity properties in neutron stars have not yet been well investi- 
gated. The aim of this letter is devoted to the influence of three-body forces on the 
1 S'o neutron and proton superfluid phases in /^-stable matter. 

For such a purpose, we shall not go beyond the BCS framework. In this case, the 
pairing gap which characterizes the superfluidity in a homogeneous Fermi system is 
determined by the standard BCS gap equation [21], i.e., 

A £ = -E^MO^A,-,, (1) 

g P 

where v(k, k') is the bare AW interaction in momentum space, E% = y^eg — ep) 2 + A|, 
ej: and ep being the s.p. energy and its value at the Fermi surface, respectively. 

In the BCS gap equation, the most important ingredients are the AW interaction 
v(k,k') and the neutron and proton s.p. energies in /3-stable matter. For the AW 
interaction, we adopt the Argonne A 1/18 two-body interaction [22] plus a microscopic 
three-body force ( TBF ) [17]. The TBF is constructed self-consistently with the A 1/18 
two-body force by using the meson-exchange current approach [18] and it contents the 
contributions from different intermediate virtual processes such as virtual nucleon- 
antinucleon pair excitations, and nucleon resonances ( for details, see Ref. [17]). The 
TBF effects on the equation of state ( EOS ) of nuclear matter and its connection to 
the relativistic effects in the DBHF approach have been reported in Ref. [18]. The 
influence of the TBF on Landau parameters in nuclear matter and on the neutrino 
mean free path in neutron stars has also been explored [23]. 

The proton fraction and the s.p. energies in /3-stable matter are calculated by 
using the BHF approach for isospin asymmetric nuclear matter [24]. In solving the 
Bethe-Goldstone equation for the G-matrix, the continuous choice [25] for the auxil- 
iary potential is adopted since it provides a much faster convergence of the hole-line 
expansion than the gap choice [26]. The effect of the TBF is included in the self- 
consistent Brueckner procedure along the same lines as in Refs. [17, 18], where an 
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effective two-body interaction is constructed to avoid the difficulty of solving the full 
three-body problem. A detailed description and justification of the method are dis- 
cussed in Refs. [17, 18]. Here we simply write down the equivalent two-body potential 
in r-space 

(nfalW'rV) = jTr £ / df^^f^r^l - v(r' 23 ))(l - v(r' 13 )) 

4 n J 

xW 3 (Ffa'r 3 '\mr 3 )(l - V(ris))(l ~ r/(r 23 )) , (2) 

where the trace is taken with respect to the spin and isospin of the third nucleon, 
and r](r) is the defect function. According to Eq. (2) the effective two-body force is 
obtained by averaging the three-body force over the wave function of the third nucleon 
taking into account the correlations between this nucleon and the two others. Due to 
its dependence on the defect function the effective two-body force is calculated self- 
consistently along with the (7-matrix and the auxiliary potential at each step of the 
iterative BHF procedure. 

With the obtained EOS of asymmetric nuclear matter the proton fraction Y p for 
a given total baryon density ps in /3-stable matter can be calculated according to the 
charge- neutrality and the equilibrium condition with respect to weak interaction [18]. 
The calculated proton fractions are listed in Tab. 1, where the first column gives the 
total baryon densities, the second column and the third column present the corre- 
sponding proton fractions obtained by using the AVI 8 two-body interaction and the 
AVA8 plus the TBF, respectively. Inclusion of the TBF in the calculations increases 
the proton fractions at high densities. This has to be attributed to the influence of 
the TBF on the EOS of asymmetric nuclear matter as verified in Ref. [18], where the 
possible implications for the neutron star cooling mechanisms were also discussed. 

To solve the gap equation, we follow the scheme given in Ref. [11], where it is shown 
that the gap equation can be split into two coupled equations, 

A k = ~ J2 V&k'^Aj:,, (3) 
4 



V(k, k>) = V(k, k>) - £ H^vfi, k>) , (4) 

k">k c ltj k" 

where the effective interaction V, arising from the introduction of a cutoff k c in mo- 
mentum space, sums up a series of ladder diagrams analogous to the Bethe-Goldstone 
equation and it is quite sensitive to the tail (A; > k c ) of the NN interaction, which 
reflects the short-range part of the nuclear force. 

Table 1 Proton fractions Y p in 5-stable matter. 
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0.0731 


0.300 


0.0667 


0.0944 


0.340 


0.0746 


0.1162 


0.400 


0.0872 


0.1499 


0.450 


0.0970 


0.1785 



In order to numerically investigate the effect of the TBF we have solved the gap 
equation, by adding the effective TBF given in Eq.(2) to the bare AV18 two-body 
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force. At the same time the s.p. energy spectrum appearing in the gap equation is 
computed from the BHF approach by using the ^4yi8 plus the same TBF. 

Fig. 1 shows the neutron energy gap in the 1 S partial wave channel A^ = A(fcp) 
as a function of the total baryon density ps- The dashed curve is obtained by adopting 
the pure two-body interaction only, while the solid curve is predicted by using 

the ^4^18 plus the TBF. For comparison, the dotted and dot-dashed curves are the 
energy gaps reported in Ref. [11] using the Argonne AV1A and the Paris potentials, 
respectively. It is seen that the results calculated with the three different two-body 
interactions agree well with each other, i.e., the neutron superfluidity phase in the 1 S 
channel can only occur in the low-density region (p# < 0.1fm~ 3 ) of neutron stars with 
a maximal gap value of about 2.8MeV peaked at a Fermi momentum k F ~ 0.8fm _1 
(the corresponding total baryon density is p# ~ 0.02fm -3 ). The TBF effect is quiet 
small, i.e., almost negligible at relatively low density and a slight suppression of the gap 
as increasing density. This result is expected from the low density for the 1 S'o neutron 
superfluidity, since three-body forces are invented to take, in an effective way, the non- 
nucleonic degrees of freedom in nuclear medium into account and become significant 
only at high densities, i.e., around and above the empirical saturation density [19]. 

In Fig. 2 is reported the proton 1 S energy gap in /3-stable matter. The solid curve is 
obtained by including the TBF while the dashed one by using the pure ylV^18 two-body 
force only. Without the TBF, our result is in good agreement with the predictions by 
adopting the BONN potential [12, 13], the AVU and the Paris potentials [11]. As 
compared to the neutron 1 S'o superfluidity, the proton 1 S superfluid phase extends to 
much higher densities but with a smaller peak gap value around pb ~ 0.2fm~ 3 . The 
former is a direct consequence of the small proton fraction in /^-stable matter, and the 
latter stems from the different s.p. potentials for neutrons and for protons [11]. As is 
known, in isospin highly asymmetric nuclear matter like /3-stable matter, the proton 
s.p. potential is much deeper than the neutron one [18, 20]. 

The effects of the TBF are twofold as shown in Fig. 2. One is a strong reduction 
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of the peak value of the gap from ~ 0.95MeV to ~ 0.55MeV and a remarkable shift 
of the peak to a much lower baryon density from ~ 0.2 fm~ 3 to ~ 0.09fm~ 3 . The 
another is that the TBF leads to a noticeable shrinking of the density region of the 
superfluid phase from p B < 0.45fm~ 3 to p B < 0.3fm~ 3 . The above predicted TBF 
suppression of the 1 S'o proton superfluidity in /3-stable matter appears inconsistent with 
the small proton fractions which correspond to small proton densities in the matter. 
However, since proton pairs are embedded in the medium of neutrons and protons, 
both the surrounding protons and neutrons contribute to the TBF renormalization of 
the proton-proton interaction. This means that the relevant density to the TBF effect 
is the total baryon density, but not the proton one. One can verify from Fig. 2 that the 
strongest suppression of the energy gap is mainly in the region ps > po, po being the 
empirical saturation density of nuclear matter, and the reduction of the gap increases 
rapidly as increasing the total baryon density. 

It is worth noticing the discrepancy between the TBF effect and the relativistic 
effect [13] on the 1 S proton superfluidity in /5-stable matter. The investigation of 
Ref. [13] shows that the relativistic effect reduces remarkably the 3 P 2 neutron gap, 
but the differences between the relativistic and nonrelativistic proton gaps in the 1 S'o 
channel are quite small. One possible reason concerns the more or less different mech- 
anisms involved. In the DBHF approach [19, 27], the medium renormalization of the 
bare NN interaction are taken into account via the Dirac spinor which is dressed in 
nuclear medium (which can be traced to the virtual excitation of nucleon-antinucleon 
pairs [28]), but the most important effect on the pairing gaps ( for instance, the 3 F 2 
neutron gap ) comes from the relativistic modification of the s.p. energies [13]. 

The TBF is expected to influence the superfluidity phases in /3-stable matter via 
three different ways. First it renormalizes the bare nucleon-nucleon interaction. Second 
it modifies the s.p. energy spectrum in the gap equation and finally the inclusion of 
the TBF changes the predicted proton fractions as shown in Tab. 1. To see which 
is the most important mechanism responsible for the strong suppression of the 1 S'o 
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proton energy gaps, we include the TBF in the BHF calculations to obtain the proton 
fractions and s.p. energies, but adopt only the pure AV18 as the NN interaction 
v(k,k') in the gap equation, i.e., Eq.(l). The results are shown by the dotted curve 
in Fig.2. One can see that the combined effects via the proton fractions and the s.p. 
energies are relatively small, i.e., a slight reduction, and mainly in the higher density 
region (pB > 0.3fm -3 ) where the TBF modifications of the proton fractions ( Tab.l 
) and the proton s.p. energies [18] become appreciably larger. Hence, the strongest 
effect stems from the TBF renormalization of the NN interaction in the medium. 

In summary, we have investigated the influence of the TBF on the neutron and 
proton pairing gaps in the : 5o channel in /3-stable neutron star matter. It is shown 
that the TBF has only a weak effect on the neutron 1 S superfluidity phase, i.e., a 
slight reduction of the energy gap, due to the low-density region concerned. However it 
suppresses strongly the proton superfluidity in the 1 S'o channel induced by the two-body 
NN interaction. The peak value of the proton 1 Sq energy gap is reduced by about 50% 
from ~ 0.95MeV to ~ 0.55MeV and shifted to a much lower density by the inclusion 
of the TBF. The density region for the superfluid phase is also remarkably shrunken 
as compared to the pure two-body force prediction. It is shown that this suppression 
is mainly related to the TBF renormalization of the 1 S'o two-body interaction. 

Besides the TBF effects, the medium renormalizations of the NN interaction and 
the s.p. energies, i.e., the screening effects such as the polarization effects and dispersive 
effects, may also influence largely the superfluidity properties in nuclear medium. Up 
to now, all investigations in the literaturefl, 2, 15, 16] have predicted a reduction of 
the BCS superfluidity gap in the 1 S channel. Therefore we expect that the screening 
effects may further suppress the 1 So proton superfluidity in /3-stable matter. 

Since the relativistic effects lead to a strong suppression of the 3 P2 neutron super- 
fluidity in /5-stable matter [13], the 1 S'0 proton superfluidity becomes critical in deter- 
mining the outcome of neutron star cooling. It is expected in Ref. [13] that the main 
suppression (of the neutrino production) comes from the superfluid proton in the 1 Sq 
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state. In this letter, the predicted suppression of the 1 S'o proton pairing gap does not 
favor the proton superfluid phase which is expected to suppress the modified URCA 
processes in the interior of a neutron star. This is compatible with the recent result of 
Link [7] derived from the observations of long-period precession in isolated pulsars. 
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Figure Captions: 

Figure 1: Neutron 1 S'o energy gap in /3-stable matter. 



Figure 2: Proton 1 S pairing gap in ^-stable matter. The solid curve is predicted by 
using the ^41^18 plus the TBF, and the dashed curve by using the Ayi8 two-body 
force only. The dotted curve is calculated with the pure Al/18 for the NN interaction 
in the gap equation and with inclusion of the TBF in the BHF calculation. 
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